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In this paper, we study the conductance and shot noise in transport through a multi-site system in 
a two terminal configuration. The dependence of the transport on the number of atoms in the atomic 
wire is investigated using a tight-binding Hamiltonian and the nonequilibrium Green's function 
method. In addition to reproducing the even-odd behavior in the transmission probability at the 
Fermi energy or the linear response conductance in the normal-atomic wire-normal metallic(NAN) 
junctions, we find the following: (i) The shot noise is larger in the even-numbered atomic wire than 
in the odd-numbered wire, (ii) The Andreev conductance displays the same even-odd parity effects 
in the normal-atomic wire-superconducting(NAS) junctions. In general, the conductance is higher 
in the odd-numbered atomic wire than in the even-numbered wire. When the number of sites (N) 
is odd and the atomic wire is mirror symmetric with respect to the center of the atomic wire, the 
conductance does not depend on the details of the hopping matrices in the atomic wire, but is solely 
determined by the coupling strength to the two leads. When N is even, the conductance is sensitive 
to the values of the hopping matrices. 

PACS numbers: 73.63.Nm, 74.50.+r 



I. INTRODUCTION 

Electron transport through an atomic wire is realized 
in scanning tunneling microscope(STM) experimentad'EI 
and break junctionsEI. An array of quantum dotscl can be 
considered as an artificial atomic wire. In the STM ex- 
periments, a true atomic wire was realized as the dipped 
gold STM tip is slowly retracted from the Au substrate. 
Just before breaking, the metallic contacts become thin- 
ner down to one atom wide and consisting of a few bridg- 
ing atoms. Yanson and othersa succeeded in forming 
chains of gold atonis—up to seven atoms long. The ob- 
served conductancea'Bcj stayed close to one unit of the 
conductance quantum before breaking. 

A quantum dot is a collection of electrons within a 
small space and can be considered as an artificial atom. 
Since electrons are confined to the small region, electrons 
experience the strong Coulomb repulsion and their en- 
ergy level spectrum is discrete. Like electrons in atoms, 
the charge and spin of electrons are quantized in quan- 
tum dots. In contrast to the atoms in the periodic table, 
it is easy to control the energy level spacing, the strength 
of Coulomb interaction, the number of electrons (charge), 
spin, etc in quantum dots. When quantum dots are ar- 
ranged as a linear chain, the system can be considered 
as an "atomic" wire. Some I-V curves have been mea- 
sured in an array of quantum dots. The splitting of the 
Coulomb peaksa was observed in the conductance when 
the interdot tunneling rate is strong compared with the 
coupling to external electrodes. 

In this paper, we study the conductance and shot noise 
in transport through an atomic wire in a two termial 
configuration. The atomic wire in our study consists of 
the monovalent atoms like Na, Au, Ag, Cu, etc. and 



is modeled by a single orbital tight-binding Hamilto- 
nian(see Fig. |l|.) When two metallic leads, connected 
to the atomic wire, are both normal, we study the de- 
pendence of linear response conductance and shot noise 
on the number of atoms in the atomic chain. When one 
lead is normal and the other is superconducting, the de- 
pendence of the Andreev conductance on the number of 
atoms in the chain is studied. Our discussion is con- 
fined to the linear response regime or to the case of weak 
source-drain bias voltage. 

The even-odd behavior of the linear response conduc- 
tance or oscillation of the conductance as a function of 
the number of atoms in the atomic wire, .was theoreti- 
cally studied using ab initio .calculations&Q. In hetero- 
geneous carbon atomic wiresu where the carbon atomic 
wires are connected to two metal electrodes, the linear 
response conductance is found to be oscillating as a func- 
tion of the number of C atoms in the atomic wire. The 
conductance is smaller than the conductance quantum 
Gq = 2e 2 jh and is higher in the odd-numbered atomic 
wires than in the even-numbered systems. On the other 
hand, the conductance is wry close to Gq in the homoge- 
neous sodium atomic wireQ when the number of atoms in 
the wire is odd. In general, the conductance is higher for 
the atomic wire with an odd number of atoms than with 
an even number of atoms. This even-odd parity effect in 
the conductance can he. understood based upon the tight- 
binding Hamiltoniarotl for the atomic wire. In addition 
to reproducing the even-odd behavior of the conductance 
in the metal-atomic wire-metal junction, we also find the 
same oscillation in the Andreev conductance as a func- 
tion of the number of atoms in the atomic wire for the 
metal-atomic wire-superconducting junctions. When the 
number of sites (N) is odd and the atomic wire is mirror 
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symmetric with respect to the center of the atomic wire, 
the conductance does not depend on the details of the 
hopping matrices in the atomic wire, but is solely deter- 
mined by the coupling strength to the two leads. When 
N is even, the conductance is sensitive to the values of 
the hopping matrices. 

In recent STM experiments of gold nanocontactsi, 
a small variation of conductance within the region of 
a unit conductance quantum was observed (Fig. la of 
referencetl) as the STM tip is stretched before breaking. 
Since the conductance is close to the unit conductance 
quantum, the gold chain is one atom thick. A small vari- 
ation but abrupt jump in the conductance might origi- 
nate from an addition of an extra gold atom to the chain. 
That is, the observed conductance variation may be a sig- 
nature of the even-odd parity effects in the atomic wire. 

This paper is organized as follows. In Sec. II, we 
study the transport through multi-sites which are con- 
nected to the two normal metallic leads. We consider 
the normal-superconducting junctions passing through 
the multi-sites in Sec. [II. Two sections are the sepa- 



Au, Ag, Cu) can be described by the model Hamiltonian, 



rate works and can be read independently. Our study is 
summarized in Sec. IV 



II. NORMAL-NORMAL METALLIC (N-N) 
JUNCTIONS CONNECTED BY MULTI-SITES 



In this section we study the current and its noise in the 
transport through multi-resonant sites in a two-terminal 
configuration where two leads are in the normal metal- 
lic state. The formalism is most conveniently simplified 
using the matrix notations. The-cucrent is given in a well- 
known Landauer-Biittiker forrrjl3lld and the.shot noise is 
also expressed in the familiar binomial formt3 except for 
some special cases. 

Applying the formalism to the serial multi-site system 
or the atomic wire connected to two metallic leads, we 
give a simple explanation of the even-odd behavior in the 
transmission probability at the Fermi energy or linear re- 
sponse conductance which was reported in the recent the- 
oretical worksoQ. In a mirror symmetric atomic wire, the 
transmission probability at the Fermi energy takes the 
same value which depends only on the coupling strength 
to the two leads, when the number of sites is odd. On the 
other hand, for the system with an even number of sites, 
the transmission probability depends on the detailed val- 
ues of the hopping matrices between neighboring sites. 
This even-odd behavior can be understood in terms of the 
level-splitting in quantum physics and the particle-hole 
symmetry in the energy spectrum when all the energy 
levels of atoms lie at the Fermi energy. We also find that 
the shot noise in the transport through multi-site system 
is oscillating as a function of the number of atoms in the 
atomic wire. 

The multi-site system made of monovalent atoms (Na, 
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Since the electron spin does not play any role in the N- 
N junctions, the spin index will be suppressed in this 
section. Electrons can flow from one lead to the other one 
through all the sites when a finite source-drain voltage is 
applied. All the sites are connected to each other via 
tunneling. The serial multi-site model is recovered when 
Wij — unless j — i ± 1 and Vn, — VzA ; i and Vm = 
VrS^n- The algebra is most conveniently simplified and 
becomes compact when the Hamiltonian is written in a 
matrix form by introducing 
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and Hd, t u = Ei,H d ,ij = Wij. With these notations, the 
model Hamiltonian becomes 
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A. Current and shot noise 

The current operator which is defined as the variation 
per unit time of electron charges in the left or right lead 
can be written in a matrix form as 



h = 



±[N p , H] = ^Y. [ c tkV^d *], W] (-2-8) 



After thermal averaging the current operator, current is 
written in terms of the mixed Green's function out of 
equilibrium. 

I p = 2e£lmTrG< (i.At)^, (2.9) 
k 

ihG dp (t,kt') = (T c y d (t)cl k (t')). (2.10) 

Here T c means contour - rOr . d ering where the contour en- 
closes the real-time axisESE-X Using the Dyson equation 
for the mixed Green's function G dp and the Green's func- 
tions of two leads, the above expression of current is i f% 
duced to a simple compact Landauer-Biittiker formt2lliil. 

I L = m deT(e)[f R (e)-f L (e)], (2.11) 



T(e) = 4ReTr[T L D r (e)T R D a (e)}. 



(2.12) 
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The factor 2 is recovered in II to take into account two 
directions of an electron spin. T p — nN p V p V^ (p = L, R) 
is the linewidth matrix. D r ' a is the retarded (advanced) 
part of Green's function of the multi-site system which 
is defined by 



iHD(t,t') = (T c * d (^(0>. 



(2.13) 



fp(p — L,R) is the Fcrmi-Dirac thermal distribution 
function in the left and right leads. 

The current noise is defined as the current-current cor- 
relation function. 

s(t,t') = (si(t)si(t')) + (si(t')si(t)). (2.14) 

Here 5I(t) = I(t) — (I(t)) is the current fluctuation op- 
erator. To compute the current noise, we introduce the 
current Green's function, 



ihG n {t,t') = (T c 5i{t)5I(t')). 



(2.15) 



Using the greater and lesser current Green's functions, 
G>(t,t') = (Si(t)SI(t')) and G<(t,t') = -(Si{t')Si(t)), 
the current noise can be written as 



S{t,t') 



hG<(t',t). 
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Since G-j Z (t,t') = — Gfj(t',t), calculation of only the 
lesser current Green's function is needed to find the cur- 
rent noise. For the noninteracting electrons in the multi- 
site system, there are only six Feynman diagrams, as 
shown in Fig. ^, contributing to the current noise. De- 
tailed calculations of each diagram are described in the 
Appendix |A|. The current noise expression is given by 
the equation, 
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Additional factor 2 is multiplied to get the spinful result 
of S . Here T = AY L D r T R D a and f p (p = L, R) is the 
Fermi-Dirac thermal function in the left and right elec- 
trodes and f p = 1 — f p . As expected, the current noise 
So is symmetrical in the lead indices L and R. In the ex- 
pression of current, the transmission probability T{e) is 
defined as T(e) = TrT. Manipulating the thermal func- 
tions, we can rewrite the expression of So as 
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The first line is the thermal or Johnson noise and re- 
mains finite at nonzero temperature. The second line is 
nonzero only out of equilibrium and is the so-called shot 
noise deriving from the discreteness of electron's charge. 



We note that the integrand in the second line is simpli- 
fied to the familiar binomial expression!! 2 ] T(e)[l — T(e)] 
except when both detr^ ^ and detr^ ^ 0. This point 
can be easily proved by noting that T and T 2 can be 
diagonalized at the same time and the trace is a sum of 
diagonal elements. Only one diagonal element is non-zero 
after diagonalization and all the others are zero. As con- 
crete examples, the familiar binomial expression of So is 
obtained for the two-site or the serial multi-site systems. 
Note that for a 2 x 2 matrix A we have the identity 



TrA 2 = [TrA] 2 - 2detA 



(2.19) 



However, detT = since detr_L = = dctr_R. That is, 
the current noise for the two-site system is reduced to 
the familiar binomial expression. There are exceptional 
casesEj of two-site system for which detr^ ^ ^= detT R . 
In this case, the familiar binomial expression of the shot 
noise is not recovered. 



B. Atomic wire: serial multi-sites 

In this section, we apply the formalism of the previ- 
ous section to the metal-atomic wire-metal junctions for 
which the atomic wire is modeled by the tight-binding 
Hamiltonian. In a serial configuration of multi-sites, the 
hopping matrices are non-zero only between two neigh- 
boring sites, 



Wij = tid i+ ij + tjSij+i. 



(2.20) 



Since only the left (right)-most site is connected to the 
left (right) lead, the linewidth matrices are simplified as 



F R ,ij = T R di^Sj^N- 
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Under this condition, the transmission probability is 
given by the expression, 

T(e) - 4TrT L D r T R D a =4T L T R |D^(e)| 2 .(2.23) 

Here D\ N is the retarded part of the Green's function 
defined by 

ihD 1N (t,t') = (T c di(*)<4 (*')>, (2.24) 

and is given by the (1, N) element of D r (e) = [e — Hd + 
tT] -1 . Here r = Tl + T R is the total linewidth matrix. 
Using the cofactor of the matrix, we find D r 1N , 



D r 1N (e) = [ e -H d + iT}^ 
tit2---tjsr-i 
~ Z N (e) ' 
Z N (e) = det[e - H d + iT]. 



(2.25) 
(2.26) 



The determinant Z^ is calculated iteratively. 

Consider the case of E t = for all i = 1,2, • • • , N 
and calculate the transmission probability at the Fermi 
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energy. This particle-hole symmetric energy structure 
seems to be realized in the ab initio calculation of the 
homogeneous sodium atomic wired. Depending on the 
parity of the number of sites, the determinant takes on 
different forms, 
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Accordingly, the transmission probability has the even- 
odd parity dependence: 
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The conditions of perfect transmission can be readily 
found from the above equations. 

When the hopping matrices are mirror symmetric with 
respect to the the center of the atomic wire, the above 
expression of Jjv(O) for odd N is further simplified. 
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Note that this value does not depend on the detailed 
values of the hopping matrices between neighboring sites. 
On the other hand, the transmission probability Tjv(0) 
for even N depends on the detailed values of hopping 
matrices. 

When the hopping matrices between neighboring sites 
are all equal, U = t for all i's, the transmission probability 
Tjv(0) is further simplified to 
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The transmission probability at the Fermi energy in the 
even-numbered atomic wire still depends on the hopping 
matrix. Perfect transmission is achieved when Tl = T R 
in the odd-numbered atomic wire, and when t 2 = r^r^ 
in the even-numbered atomic wire. 

The even-odd parity behavior of the transmission prob- 
ability at the Fermi energy is a consequence of the energy 
level splitting in quantum physics and the particle-hole 
symmetry of atomic energy levels. In the particle-hole 



symmetric case, all the atomic energy levels lie at the 
Fermi energy (Ei = 0) when the hopping between atoms 
and the coupling to the leads are turned off. When hop- 
ping between neighboring atoms in the atomic wire is 
turned on, the degeneracy in atomic energy levels is lifted 
and the energy level structure still keeps the particle-hole 
symmetry (see the Appendix |b|.) Note that the values of 
tj's have no influence on the particle-hole symmetry of 
split atomic energy levels. 

In the odd-numbered atomic wire, one level always lies 
at the Fermi energy and all others are distributed sym- 
metrically with respect to the Fermi eneigy(Ep). Note 
that Z 2 n+i(0) is equal to when Tl, R = 0. This vanish- 
ing determinant means that at least one energy level lies 
at Ep. Since no degeneracy remains in the presence of 
hopping between atomic sites, only one level is located at 
the Fermi energy. When the coupling to external leads is 
turned on, the discrete energy levels become broadened. 
The transmission probability T/v(e) for odd N is always 
peaked at e = or the Fermi energy. 

On the other hand, all energy levels are distributed 
symmetrically with respect to the Fermi energy without 
any level at Ep for the even- numbered atomic wires. 
Since Z 2 n{0) is not equal to when = 0, the 

degeneracy-lifted atomic energy levels cannot lie at the 
Fermi energy in the even-numbered atomic wire. Though 
no energy level is present at the Fermi energy for the 
even-numbered atomic wire, either a dip or a peak in 
Tjv(e) for even N is possible at the Fermi energy e = 
when Tl R ^ 0, depending on the relative magnitudes of 
the hopping integral (level-splitting) in the atomic wire 
and the coupling strength (linewidth) to the leads. When 
ti = t for all i's and Tl = T R , the transmission probabil- 
ity is dipped (peaked) at the Fermi energy(e = 0) when 
t > {<)Tl,r, respectively. 

To illustrate the even-odd parity dependence of the 
transmission probability, Tjv(e) is displayed in Fig. || for 
differing N's (the number of atoms in the atomic wir. 
To reproduce the results of the ab initio calculations 
model parameters are chosen such that Ei = for all 
sites and all tight-binding hopping matrices are equal, 
ti = t. For typical metals, t is of the order of leV. Fur- 
thermore, the bonding in the gold jatomic wire is almost 
two times stronger than in the bulkB. To simulate the ho- 
mogeneous monovalent atomic wires (Au, Na, Ag, Cu), 
we take some intermediate value oiTp ji/t = 0.7 in Fig. |3|. 
It is a reasonable approximation to assume that the bond- 
ing between the atom in the atomic wire and the atom 
in the leads is intermediate compared to the bonding in 
the chain or in the bulk. Clearly, T(e) is peaked at the 
Fermi energy in the odd-numbered atomic wires and is 
suppressed at e = in the even- numbered atomic wires. 
The number of peaks in T(e) agrees with the number of 
atoms in the atomic wire. 

Sharp asperity at the STP tip seems to be necessary 
for the conductance quantization according to recent nu- 
merical workal3. In our approach, the blunt STM tip 
can be simulated by increasing the coupling strength be- 
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tween the electrodes and the atomic wire. Note that the 
left (right)-most atom in the chain is coupled to more 
atoms in the leads for blunt tips than for sharp tips. The 
results for TL,B,/t = 1.2 arc displayed in Fig. [|. Though 
the spectral shape of 7jy(e) is broadened compared to 
Fig. |3| the even-odd parity effects persist in this case, 
too. Note that the dip of T2jv(e) at the Fermi energy for 
t > ^l.r is transformed into a peak when t < Cl._r. 

In the heterogeneous carbon atomic wirescl, charge 
transfer occurs from the metallic electrode into the car- 
bon chain, and a Schottky-like barrier is formed be- 
tween the electrodes and the carbon chain. Due to the 
Schottky-like barrier, the coupling strength between the 
electrodes and the chain is weaker than in the homoge- 
neous systems. The authors suggest that the reduced 
transmission probability, T(0) < 1, is a result of the 
higher barrier between the chain and the leads. Higher 
barrier or weaker coupling leads to sharper line shape in 
the transmission probability. We would like to propose 
a different interpretation. As long as the couplings to 
two electrodes are symmetrical and all the hopping ma- 
trices are the same, the broken particle-hole symmetry 
in the energy levels of atoms in the chain can be the rea- 
son for the reduced T(0) < 1. The doped excess charge 
is not distributed uniformly over C atoms in the chainH. 
The assumption of Ei = for all atomic sites seems to 
be violated. According to the Friedel sum rulcO, ex- 
cess (depleted) charges will shift downward (upward) the 
energy levels B, of atoms in the chain. Broken particle- 
hole symmetry in the energy structure leads to reduced 
transmission probability at the Fermi energy in the tight- 
binding Hamiltonian approach. 

We can also study the effects of disorder in the hopping 
matrices in the atomic wire on the transmission proba- 
bility. While retracting the STM tip from the sample, 
atomic rearrangements occur according to molecular dy- 
namics simulationsa. Some randomness in the bonding 
between neighboring atoms might be present. To simu- 
late this situation, we introduce the disorder in the mag- 
nitude of i's such that i< = t(l + <5;) with some random- 
ness Si. Typical results of T(e) with disordered V s with 
\5i\ < 0.2 are displayed in Fig. ||. Even in the presence 
of disorder in the hopping matrices, the even-odd parity 
behavior may persist as far as the disorder is not strong 
and the coupling to the leads is not close to the values of 
hopping matrices. 

We predict that the shot noise also displays the even- 
odd behavior depending on the number of atoms in the 
atomic wire. First note that TrT 2 = [T(e)] 2 due to the 
simple structure of the linewidth matrices. Hence the 
shot noise is given by the familiar binomial expressiontL 2 ! 
and at zero temperature by 

= 4e2 p deT(e)[1 _ T(e)] _ (2 36) 

When the source-drain bias voltage is small such that 
the transmission probability does not change much over 
the energy eV = fJ-L — [Ir near the Fermi energy, the 



shot noise can be approximated as Sq = (4e 2 //i)T(Q)[l — 
T(0)]eV. In the homogeneous sodium atomic wiraJ, the 
transmission probability was found to be very close to 
unity when the number of atoms is odd. On the other 
hand, the transmission probability is smaller than unity 
when the number of atoms is even. We conclude that the 
shot noise at low source-drain bias voltages is larger in the 
even-numbered atomic wires than in the odd-numbered 
atomic wires. 



III. NORMAL-SUPERCONDUCTING 
JUNCTIONS CONNECTED BY MULTI-SITES 

In this section, we study the I-V curves in the normal- 
superconducting(N-S) junctions connected by multi-sites 
in between. The left lead is a normal metal and the right 
lead is assumed to be the simple BCS s-wave supercon- 
ductors with the constant energy gap. In our simple ap- 
proach, the energy gap in the superconducting lead is 
assumed to be constant and not degraded near the junc- 
tion by the atomic wire. Local density of states (DOS) 
at the atomic sites is modified by the superconducting 
energy gap of the right lead and its structure depends on 
the parity in the number of atoms in the atomic wire. 
In addition, the DOS structure takes different shapes de- 
pending on the atom's position in the chain. As in the 
N-N junctions connected by an atomic wire, the Andreev 
conductance shows the even-odd parity behavior. 

To describe the superconducting state, it is convenient 
to use the Nambu spinor notation. The two leads are 
described by the Hamiltonian, 

H <* = E T,Kk E p^ P ^ (3-i) 

p=L,R k 

** = (T' ) ■ (3-3) 

\ C p-klJ 

Here e p k is the energy dispersion in the left (p — L) and 
right [p = R) leads and A r is the energy gap in the right 
superconducting lead. Since the phase of the energy gap 
An in the N-S junction does not play any role, An is 
taken to be real. The Hamiltonian of the multi-sites is 
also written using the Nambu spinor notation as 

+Y,*U w o -S*)**' (3 - 4) 

where = (d\^ dijj . Ei is the energy level of the z-th 

site and Wij is the hopping matrix between two sites i 
and j. The coupling between two leads and the multi- 
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sites is described by the tunneling Hamiltonian, 



pki 



(3.5) 



Using this-, model Hamiltonian, we study the Andreev 
reflectionta in the N-S junction. 



A. Multi site: General formalism 

The algebra is highly simplified when the Hamiltonian 
of the multi-site system is written in terms of the matrix. 
Introducing new notations 
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the model Hamiltonian can be written as 
H d = tf d H d i& d , 



(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 



Here t 3 is the third component of the Pauli matrices. The 
current operator can be also written in a simple compact 
form using the new notation, 



ih 

k 



H.c. 



(3.11) 



The thermally averaged current is then 



I p = 2e^ImTrVj G<(t,kt). 



(3.12) 



The mixed Green's function is defined and determined 
by the Dyson equation, 



1 



G dp (t,kt') = ^{TcMtHkit')) 



dti D(t, h)V p T 3 G cp (k; h, t')(3.13) 



ihD(t,t') - (T c H> d {t)H> d (t')) 



(3.14) 



Inserting all the relevant Green's functions of two leads, 
the current can be expressed as 
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This equation of current is our starting point for the 
study of the even-odd parity behavior in the Andreev 
conductance in the metal-atomic wire-superconducting 
junctions. D r ' a is the retarded (advanced) part of the 
Green's function of the multi-site system. Note that 
D a = [_ D r]t_ other notations are listed below. 



Me) 
9(e) 



f(e-»L) 
f(e + »L, 





)' 


I 








7T (e-f 








-N) 






-e 2 



isgn(e) 



A R \ 2 



A, 



A 2 



(3.16) 
(3.17) 

3.18) 



The second line in g is obtained in a wide conduction 
band limit. We write g = gi + ig 2 - Note that g is real 
for |e| < A R and imaginary for |e| > A R . To simplify 
the algebra, the bias voltages are chosen as = eV and 



A'-R — 0. The first line inpEq. (3.15) is a contribution 
from the Andreev reflectionES. The integrand in the first 
line remains finite over the entire energy range, though 
it is appreciable inside the gap. The electron incident 
from the normal metallic lead forms a Cooper pair with 
its partner, tunnels into the superconducting lead and 
leave a hole in the normal metallic lead. This remaining 
hole moves along the time-reversed track of the incident 

electron. Note that .(72( e) for |e| < A^. That is, 

the second line in Eq. ( 3.15| ) is the contribution from 
quasiparticles excited over the superconducting energy 
gap. The above expression of current can be applied to 
any configuration of multiple resonant sites connected to 
the two leads. 



B. One site 

To begin we consider the N-S junction connected by 
the one-site system. The coupling matrices to the leads 
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V. 



v L1 
o v 1L 



v R 



V m 
Vi 



. (3.19) 



Here Z\ is the determinant of [-D r (e)] 1 and the denom- 
inator T> is given by the expressions, 



For the one-site system, the Eq. ( |3.15 ) becomes simplified 
to 



h 



4e, 



deTr 



D r f L D a r 3 - D r r 3 f L D a 



V = 



(e 2 - E( - T 2 L - T R )J A 



-4r?e 2 



4e f 

r D r r 3 f L D a - D a T 3 f(e)D r ] . (3.20) for lei < A R and 



R 

n 2 



r R + JA 2 -e 2 



2T R e 2 



(3.25) 



Here T p (p = L, R) is the linewidth parameter defined 
by the relation, T p — irN p \V p i\ 2 . The first line is the 
contribution from the Andreev reflection and the second 
line is the contribution from quasiparticles excited over 
the superconducting energy gap. The self-energy of a 
resonant site is given by the equation, 

£ r = -zT L l + T R g(e)T 3 Cl R T 3 , (3.21) 
S a = iT L l + T R g*(e)T 3 Cl R T 3 . (3.22) 

Here 1 is a 2 x 2 unit matrix. Inserting the Green's 
functions of a resonant site, the current can be written 
in terms of the transmission probability T(e), 



V = 



(e 2 -E 2 1 -T 2 L -T 2 R ) x /e 2 -A 2 R -2T L T R \e\ 

2 



-Ae 2 



T L Je 2 -A 2 R + T R \e\ 



(3.26) 



for |e| > A R . Though the detailed forms of T A {e) are 
different depending on the energy range, the Andreev 
contribution is continuous at |e| = A R . 



It 



2e 



deT(e) [/(e) - /(e - (x L )] . (3.23) 



T(e) is defined as a sum of two, T(e) = T A (e) + T qp (e) 
The Andreev reflection gives 



T A (e) = 8T 2 L \D[ 2 (e) 



A 2 



IZ.ieW \e 2 - A 2 R \ 



V(e) 



(3.24) 



Ta(0) 



T a (±A r ) 



{E 2 l+ T 2 L +T R y 

2 r| 



(3.27) 



(3.28) 



Note that Ta at e = ±A R is independent of the energy 
level position i?i and is always less than a value of 2. The 
quasiparticle contribution is given in an explicit form as 



T qp (e) = 4T L T R \g 2 (e)\ [|e| {^(e)! 2 + \D r l2 {e)\ 2 } - A fl sgn(e)Re {^(e^e) - ^ 2 (-e)£>? 2 (-e)}] 



4r L r i? . v /e 2 - a^ 


|e| { (e + e x ) 2 + r| + r|} + 2r L r flv / e 2 - K 




( £ 2 - e\ r 2 - i* ) Ve 2 - A 2 fl - 2r L r fl | e |_ 


2 + 4e 2 


r^-A^ + rfliei 


2 



As expected, the quasiparticle contribution vanishes at 
|e| = A R . When the right lead becomes normal or 
A R = 0, the Andreev contribution vanishes and the 
quasiparticle contribution is simplified to the well-known 
Lorentzian form, T qp = 4T L T R /[(e - E^ 2 + (T L +T R ) 2 }. 



Atomic wire: serial multi-sites 



When multi-sites are arranged in series like atomic 
wires, the current can be simplified considerably. Since 
the left (right)-most site is coupled to the left (right) lead, 
respectively, the coupling matrices to the leads are given 
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by the expressions, 



2x2 matrix in the Nambu spinor space. 



V 



v R = 



V L1 
o V 1L 



- o 



v RN 

V, 



NR 



(3.30) 
(3.31) 



Here is the 2x2 zero matrix. Various matrices appear- 
ing in the expression of current are also simplified. For 
example, 



[ttNlVlV^ = TlIS^iSj,!. 



(3.32) 



Here 1 is the 2x2 unit matrix in the Nambu spinor space. 
The Andreev contribution to current can be written in 
components as 



Ia 



I A = ^ jdeT A (e)[f(e)-f(e-f, L )}, (3.34) 
T A (e) = 4r|[|^ M2 ( £ )| 2 + |^ 142 (- e )| 2 ].(3.35) 

The Andreev conductance is given in a simple form, 

Gns = 2-^-4ri|^ M2 (0)| 2 . (3.36) 

The quasiparticle contribution to current is also given in 
components by the equation, 



?/<* 

^ • T| J de\D r lhl2 (e)\ 2 [f(e + » L ) f(e ^3.33) 



iqp 



4e r 
~~hJ 



de r L r fl (e)sgn(e)TrT 3 O fi T3 



x W N J L r 3 T> a 1N - T> N1 fr 3 W 1N . (3.37) 



The indices i and j in Dij iCl p denote the sites of resonant 
levels while the indices a and (3 mean the elements of the 



Here r#(e) = | <7 2 C e ) I - After some algebra, the quasi- 
particle contribution to the current can be written as 



T qP (e) 



I an - 



| J deT qp (e) [/(e) - f(e » L )] , 

2T L T R (e)\e\ [\D Nhll (e)\ 2 + |D^, 12 (e)| 2 + \D Nh22 (-e)\ 2 + \D N1A2 (-e)f 
-4r L r K (e)A«sgn(6) Re [D r NliU (e)Dl Nil2 (e) - D Nh22 (-e)D^ 21 (-e)] . 



(3.38) 



(3.39) 



The total transmission probability is a sum of two con- 
tributions: T = T A + T qp . Though we are using the word 
"probability", it is a misnomer because T can become 
larger than a unity due to the Andreev reflection below 
the superconducting energy gap. 

In a serial multi-site configuration, the retarded self- 
energy is given by the expression, 

Sy- = -iT L 18iA?>j,i + F R g(e)T 3 £l R T 3 5 liN 5 J ^3A0) 

and the retarded Green's function of multi-site system 
is expressed by the matrix inversion, D r = [e — Hd — 
£ r ] _1 . We now calculate the Andreev conductance of 
multi-site atomic wire one by one. For one resonant site, 
the Green's function is given by the equation, 

n r (e) _ feZ R - Ex+iTL T R A R g{e) \ . . 
U ( ) ~ I r R A R g(e) eZ R + E 1+ iT L ) 



Gns = 2 



4T 2 L T 2 R 



2e 2 



h [Ef + Tl + rlY 



(3.42) 



Here Z R = 1 — T R g(e) and the one-site system is detailed 
in the previous section. For two-site atomic wire, the 



Green's function relevant to the Andreev reflection is 

\h\ 2 A R T R g(e) 



Z2(e) 



(3.43) 



where 



Z 2 (e) = [(e + iT L ) 2 -E 2 }[e 2 -E 2 -T 2 R -2e 2 T R g] 
-2\h\ 2 [(e + iT L )(e - eT R g) + E X E 2 \ 
+ |ti| 4 . (3.44) 



The Andreev conductance is 



2e 2 4r 2 r 2 ^ 1 ' 4 



Gats = 2 • —r- ■ 



(3.45) 



h |^i(0)| 2 ' 
Z 2 (0) = (E 2 + T 2 L )(E 2 2+ T 2 R ) + \t^ 

-2\t 1 \ 2 E 1 E 2 . (3.46) 

At e = ±A R , T A is independent of E 2 and is given by 
the expression, 



T A = 



2T\t\ 



A 2 R (E 2 +r 2 L +t 2 -A 2 R r+T 2 L (2A 2 R -4) 



#47) 
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Ta at e = ±A.r can reach the maximum value of 2 when 
Ei = and = \/t 2 — A 2 ^. For three-site atomic wire, 
the Green's function at the Fermi energy is 



where 



^11,12(0) 



\h\ 2 \t 2 \ 2 T R 



(3.48) 



case of the mirror symmetric chain), the transmission 
probability at the Fermi energy is simplified as 



T 2J v+i(0) 
T 2N {0) 



8-r| 



r 2 

1 R 



n 



8 • v L v R • t 



(3.58) 
(3.59) 



Z 3 (0) 



det(-H d + iT) 

-{El + r 2 ^ 2 
-M 4 (£ 3 2 



l) E 2 ( E 3 



T%)-\t 2 \HEf+Tl) 
2E 2 E 3 \t 2 \ 2 (E 2 + T 2 L ) 



+ \h\ 2 {E 1 E 2 {2E 2 + T| ) - £i£ 3 |i 2 | 2 } 
-i\h\ 2 T L (E 2 T 2 L -E 3 \t 2 \ 2 ). (3.49) 

When Ei = for all i's, the Andreev conductance is given 
in a simple form, 



£11,12(0) 



M 2 N 2 r 



R 



r|M 4 

2e 2 

Gns = 2 • — — • 



4r||t 2 | 4 .r| i |i 1 | 4 



n\t 2 \ 4 



r|M 4 ] 



(3.50) 
(3.51) 



We can now deduce the general expression of the Andreev 
conductance for TV-site atomic wire. Depending on the 
parity of the number of atoms in the atomic wire, the 
Andreev conductance takes different forms, 



Gns 
71(0) 

7W(0) 

T 2 (0) 

T 2N (0) 



2e 2 

^-•7V(0), 



[r 



r 2 1 

1 R\ 



(3.52) 
(3.53) 



•r 2 k 2 t 4 - 



• t 2 N\ A ■ v R \t\t^ • • • t 



2N- 



[T 2 L \t 2 U---t 2N \ i + T 2 R \t 1 t 3 - 

8-r|r 2 - 1*! | 4 



■ t 2 N-: 



[r|r 2 + M 4 ] 2: 



3^ 



(3.55) 



L 1 - R 

■Tl\t 2 U- 



■ t 2 N- 2 \ i ^ 2 n 



[T 2 L \t 2 t i ---t 2N ^T 2 R + \t 1 t 3 



tit 3 ---t 2N ^ 



■ -t 2 N-l\ 



The even-odd behavior is also evident in the Andreev 
conductance as in the NAN junction connected by the 
atomic wire. Since the Andreev reflection is the two- 
particle process, all the parameters in the transmission 
probability at Ep in the NAS junction are squared com- 
pared to the NAN junction. In the mirror symmetric 
atomic wire, T 2 n + i(0) does not depend on the detailed 
values of the hopping matrices in the atomic chain but is 
determined by the couplings to the leads. 



T 2JV +i(0) 



r 2 

1 R 



[r 



(3.57) 



T^aKO) depends on the hopping matrices in the atomic 
wire. When all the hopping integrals are equal (a special 



The maximum transmission is achieved when Tl = Tr 
in the odd-numbered atomic wire, and when t~ = TiXr. 
in the even-numbered atomic wire. 

To illustrate the even-odd parity dependence of the 
transmission probability T/v(e) = Ea{^) +T qp (e), we plot 
T;v(e) for differing TV's (see Figs. | and |). The model 
parameters are chosen as follows: Ei — for all the sites 
and ti — t for all z's. The superconducting energy gap 
Ar = t/10 is chosen to be large to show clearly the 
Andreev structure in Tjv(e). The results for the case 
Fl.r = 0.3 x t are displayed in Fig. ^. Note that the An- 
dreev contribution of Ta (e) is appreciable only inside the 
superconducting energy gap. For odd-numbered atomic 
wires, the transmission probability displays the Andreev 
peak at the Fermi energy and the side peaks above the 
superconducting energy gap which are the quasiparti- 
cle contributions. For even-numbered atomic wires, the 
transmission probability is suppressed at the Fermi en- 
ergy and is peaked at the superconducting energy gap. 
Other peaks outside the superconducting energy gap are 
contributions from the quasiparticles. 

When the coupling to the leads is larger than the hop- 
ping integral in the atomic wire, the spectral shape of 
T/v(e) is different compared to the case of < t. 

The results of Tjv(e) when T^ r — 1.2 x t are displayed 
in Fig. 0. Increased coupling to the leads results in a 
broadened spectral shape of the transmission probabil- 
ity. In the odd-numbered atomic wire, the transmission 
probability is almost flat inside the superconducting gap. 
In the even- numbered atomic wire, the Andreev con- 
ductance is increased as is increased(£ > VTYXr), 
reaches the maximum value of 2 when t = VTlTr and 
is decreased with further increasing TL^iV^L^R > t). 

The local DOS at each atomic site in the atomic wire 
is computed and displayed in Figs. || and [| The shape 
of local DOS is dependent on the parity of the number of 
atoms in the atomic wire. The local DOS also depends 
on the atom's position in the chain. 



IV. 



SUMMARY AND CONCLUSION 



Using the tight-binding Hamiltonian, we studied the 
dependence of the conductance and shot noise on the 
number of atoms in the atomic wire when the atomic 
wire is connected to the two metallic leads or to the 
metal and superconducting leads. In metal-atomic wire- 
metal(NAN) junctions, the even-odd parity dependence 
of the conductance can be understood by the energy- 
level splitting in quantum physics and the particle-hole 



10 



symmetry of atomic energy levels. The linear response 
conductance is larger in the odd-numbered atomic wire 
than in the even-numbered atomic wire. The conduc- 
tance can reach the maximum possible value 2e 2 /h in 
the odd-numbered atomic wire when all the atomic en- 
ergy levels lie at the Fermi energy (resonance condition) 
and the system is mirror symmetric with respect to the 
center of the wire. The shot noise also displays the even- 
odd parity effects depending on the number of atoms in 
the atomic wire. In contrast to the conductance, the shot 
noise is larger in the even-numbered atomic wire than in 
the odd-numbered atomic wire. 

We also studied the Andreev conductance in the metal- 
atomic wire-superconducting junctions. As in the NAN 
junctions, the Andreev conductance displays the even- 
odd parity effects. In the odd-numbered atomic wire, the 
conductance is enhanced and close to 4e 2 /h due to the 
Andreev reflection. The maximum Andreev conductance 
4e 2 //i is possible in the mirror-symmetric odd-numbered 
atomic wire when the resonance condition is satisfied. 
On the other hand, the conductance is suppressed below 



Ae 2 /h in the even-numbered atomic wire even when the 
resonance condition is satisfied for the mirror-symmetric 
case. 
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APPENDIX A: DERIVATION OF EQ. (2.17) 



In this Appendix, we find the expression of the cur- 
rent noise by calculating the Feynman diagrams shown 
in Fig. |[ We calculate each diagram step by step. The 
diagram in Fig. ||(a) gives 



iKG u {tJ) = J c dt i J c dt * ^{D(tM)V L G cL {k'-M,t')vlD(t\t 2 )V L G cL {k-M,t)vlY (Al) 

After some algebra, the lesser part is given by the equation 

hG<(e) = & - Jd(Tr{[ 1 D<(e + + 2f L (e + OD r (e + t)]r L [iD>(() + 2f L (()D r (t)]r L }. (A2) 

We use the analytic continuation to the real-time axisQ to find the lesser current Green's function. The diagram in 
Fig. |(b) gives 

ihG n (t,t') = --^ V / dh [ dt 2 TiivLG^t^vlDit^t^VLGcLik'-t' M)vlD{t 2l t)\ , (A3) 
V kk' Jc Jc L 

HG<(e) = & - J d( Tr {T L [-iD<(e + () + 2f L (e + ()D a (e + 0} T L [-W>(0 + 2f L (()D a (()} } . (A4) 
Two diagrams in Fig. |^(c) give 

2 

iW n (t,t') = yJ2 Tr { v L G cL(k;t,t')vlD(t\t)} 

k 

+ ^ J2j c dtl J c dt2 Tr { V LGc L (k; t, t 2 )v[D(t 2 , t Y )V L G cL {k'; t u t')VlD{t' , t)} , (A5) 

hGfj(e) = J dC/4e + C)Tr{r L £>(C)} 

~ I d(Tr{r L [D<(e + () + 2 t .f L (e + 0{D a (e + 0-D r (e + 0}}r L D>(()}. (A6) 
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Two diagrams in Fig. ||(d) give 



ihG n {t,t') = ^Y,^{ D ^ t ') y L G cL{k ] t\t)vl] 

k 

+ ^ J2j c dtl J c dt2 Tr { D & ^LGcUk'; t', t 2 )vlD(t 2 ,h)V L G cL {k; h,t)Vl} , 



HG<(e) 



„2 



_£_ | dC Tr {£<( e + QT L [D>(0 + 2if L (Q{D a (Q ~ D r (()}] T L } 
Collecting all contributions, the lesser current Green's function is given by the expression, 



hG<(e) 



dCTr[2/ L (e + ()D r (e + () + iD<(e + 0}T L [2f L (C)D r (0 + iD>(0]T L 

d(TT[2f L (e + ()D a (e + Q- iD<(e + Q]T L [2f L (C)D a (0 - iD>(Q}r L 
/ d C/L(e + C)Tr{r LJ D>(C)} 
~ J d( Tr {T L [D<(e + () + 2if L (e + (){D a (e + - D r (e + {)}] r L D>(()} 
~ J dC f L (()Tr{D<(e + ()r L } 

-j JdC Tr {D<(e + ()T L [D>(() + 2if L (C){D a (0 D r (C)}] T L } . 



(A7) 



(A8) 



(A9) 



The current noise is then given by 5(e) = —HGf I (e) — 
hGfj(-e). In particular, the u> = component current 
noise So = —2hGf I (0) is given by the equation, 



8e 2 



S Q = ^- / de Trr L z? < r L ^ > 



8e 



~Y J de f L fLTiT L [D r T L D r + D a T L D a ] 
~ I de f L TiT L D>T L ■ i[D r - D a ] 
~ [ def L TYT L D<T L -i[D r 



h 



D a ] 



de [f L TrT L D> + f L TrT L D<] . (A10) 



The above expression of Sq looks highly asymmetrical in 
the indices of left and right leads or L and R. Substitut- 
ing the following identity relations 



D<> = D r "E < ' > D a , 

= VlTl + 2f R T R , 

= V L T L + 2f R T R , 

i[D r -D a ] = 2D r TD a , 



(All) 
(A12) 
(A13) 
(A14) 



we simplify the expression of current noise Sq as 

s ° = T" / de[hf ~ L + /il ^ TrT2 



2e^ 
" h 



J de[f L f R + f L f R ] [TrT - TrT 2 ](A15) 
Here T = AY L D r T R D a . 



APPENDIX B: PROPERTIE S OF 
DETERMINANT Z N (e): EQ. ( |2.26D 

In this Appendix, we study the properties of th e de- 
terminant Zjv(e) which is defined by the Eq. ( 2.26| ). We 
can readily find the iterative relation for the determinant 
when Tl R — and = for all i's, 

Z N +2{e) = eZ N+1 {e)~\t N+1 \ 2 Z N {e), (Bl) 
Z 2 (e) = e 2 -\ tl \ 2 , (B2) 
Z 3 (e) = e(e 2 -|t 1 | 2 -|t 2 | 2 ). (B3) 

Especially when U = t for all i's, the expression of Zj^(e) 
can be found in an explicit form, 



Z N (e) 



N 

n 

k=l 



e — 2\t\ cos ■ 



kit 
N+l 



(B4) 
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Zpf(e) is the type-II Chebyshev polynomial. 

From the recursive relation, we can deduce that Zpj(e) 
is invariant under the transformation, ij — > — for 
any i. We can also show from the recursive relation 
that Zn(c) satisfies the parity relation: Zn{—e) = 
(— 1) ,Zjv(e). This parity relation has the important con- 
sequence on the particle-hole symmetry of split atomic 
energy level structure. Z 2 n(^) is the even function of 
e and Z 2N (0) = (-l) N \ht 3 - ■ -t 2 N-i\ 2 - Z 2N+1 {e) can 
be written as Z 2 ^ + i{e) = e^2jv+i(e) where Y 2 jy +1 (e) is 
the even function of e and Y2jv+i(0) 7^ 0. When all the 
atomic energy levels lie at the Fermi energy, the hopping 
integrals between neighboring sites in the atomic wire 
lift the degeneracy but does not change the particle-hole 



symmetry of split energy levels. If eo is a solution of 
Zn{z) = 0, so is — eo because of the parity relation. Since 
Zn(0) 7^ for even N, no split levels lie at the Fermi en- 
ergy. Since Zjv(0) = for odd N, e = is a solution 
of the secular equation Zjv(e) = 0. Since Yzv(0) 7^ for 
odd N, only one level lies at the Fermi energy for the 
odd-numbered atomic wire. 

Though the particle-hole symmetry of split atomic en- 
ergy levels is shown explicitly, it can be also proved by 
noting that the tight-binding Hamiltonian is invariant 
under the particle-hole transformation: Cj — > (— l)*c|. 

Obviously the Hamiltonian H = J2i(ti c l+i c i + H.c.) re- 
mains invariant under this particle-hole transformation. 
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FIG. 2: Feynman diagrams for the current noise. Solid line 
represents the conduction electron Green's function, while the 
dashed line means the fully dressed Green's function of the 
multi-site system. 
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FIG. 3: Dependence of the transmission probability on the 
number of atoms in atomic wires in the metal-atomic wire- 
metal junctions. Panel (a) displays the transmission proba- 
bility for the odd-numbered atomic wires. T(e) is peaked at 
the Fermi energy. On the other hand, panel (b) shows the 
transmission probability for the even-numbered atomic wires. 
T(e) is suppressed at the Fermi energy. We take the model 
parameters: Ei = and U = t for all sites, Fl,r = 0.7 x t. 
See the text for explanations. 
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FIG. 4: The transmission probability for blunt tips. Model 
parameters are the same as in Fig. |^ except for Fl,r = 1.2 x t. 
Increased coupling to the leads broadens the lineshape of the 
transmission probability. T2jv(e) is peaked at e = when 
Tl,r >t. 




FIG. 5: The transmission probability for disordered hopping 
matrices. Model parameters are the same as in Fig. ^ except 
that some disorder in ti's is introduced: i* = t(l + Si) with 
|^| < 0.2. 




FIG. 6: Dependence of the transmission probability on the 
number of atoms in an atomic wire in the metal-atomic wire- 
superconducting junctions. Panel (a) for the odd-numbered 
atomic wires and panel (b) for the even-numbered atomic 
wires. We take model parameters: Ei = and ti = t for 
all sites, Yl,r = 0.3 x t and the superconducting energy gap 
A R = 0.1 x t. 




FIG. 7: Dependence of the transmission probability on the 
number of atoms in an atomic wire in the metal-atomic wire- 
superconducting junctions. Panel (a) for the odd-numbered 
atomic wires and panel (b) for the even-numbered atomic 
wires. We take model parameters: Ei = and ti = t for 
all sites, Yl,r = 1.2x4 and the superconducting energy gap 
A R = 0.1 x t. 
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FIG. 8: Local DOS at each atomic site for 4-site atomic wire. 
Ai means the spectral function of the i-th atoms from the left 
lead (normal metals) . Model parameters are the same as in 
Fig. | 



